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1. Introduction 



A new approach to the spin \ quantum Heisenberg model with Hamiltonian: 

H= £ 4(x,y) (1.1) 

<a;,{/> 

where: 

(*,y) = i(l-*(aO*(y)) (1-2) 

has been proposed by one of us in the preceding paper [1], which is refered as paper 1 
in the following. In this second paper we are going to exploit the new method in the 
sector with total spin 0. It was pointed out in paper 1, that the partition function in 
the spin sector: 

tr (exp/3if ) = -£/(/?, (1.3) 
ao K 

contains already the whole information on the thermodynamics of the antiferromag- 
netic Heisenberg model in the absence of an external field - i.e. for vanishing mag- 
netization M = y = 0. Moreover, it was shown that the partition function (1.3) 
is obtained from the expectation values of the Boltzmann factor exp((3H ) between 
valence bond states \K): 

f(P,K) = {K\exp(PH )\K). (1.4) 
These expectation values were proven to obey an evolution equation : 

J-/(/3,X) = (dV + 3N 1 (K))f(^,K) 



dp 



+ Y, W,Q + (x,y)K)-f((3,Q_(x,y)K)) (1.5) 

<x,y> 



where Q±(x, y) are permutation operators which act on the valence bond configuration 
K, as defined in eqs.(2.13,14) of paper 1. Ni(K) is the number of dimers. These are 
the nearest neighbour valence bonds on the configuration K. 

It is clear from the evolution equation (1.5) that the expectation values f(/3,K) in- 
crease much stronger with (3 on configurations with a large dimer density x = On 
the other hand, these configurations have low entropy, since the number of valence 
bond configurations v(Ni,V) with a given number of dimers N 1 decrease as (iVx!) -1 . 
Therefore, if we assume that the averages: 

/(AiV 1 ) = — ^— £ mKiN,)) (1.6) 

V) K(Ni) 

over all valence bond configurations K(Ni) with a fixed number Ni of dimers behave 
in the combined limit: 

V,JVi,n->oo, x = N 1 /V, (3 fixed (1.7) 



as: 



f((3,N 1 )=exp(V<P(P,x)) (1.8) 
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one is led to the conclusion that the zero density (x = 0) contribution of the dimers 
determines the thermodynamical properties of the model. This statement holds in 
general provided that the function <f>(/3, x) is differentiable with respect to x. 
E.g. the internal energy per site at fixed inverse temperature (3 is given by: 

U((3) = ^((3,x = 0). (1.9) 

We are not yet in the position to evaluate the exact evolution equation (1.5) in the 
thermodynamical limit. However it will be shown in section 2 that the lowest order 
in the 'dimer cluster expansion' can be solved analytically in the thermodynamical 
limit. Results for the ground-state energy and the specific heat of the AFH- model in 
dimensions d = 1,2,3 are presented in sections 3 and 4. 



2. The Dimer Cluster Expansion in the Spin Sector 

In paper 1 an approximation scheme was proposed for the evaluation of the evolution 
equation (1.5) in the thermodynamical limit. In lowest order of this 'dimer cluster 
expansion' (1.5) is averaged over all valence bond configurations K with a fixed number 
Ni(K) of dimers. The approximation would be exact if the right hand side of (1.5) 
could be expressed as well in terms of the averages (1.6). This is not possible for all 
the terms in the sum on the right hand side of (1.5). Here we approximate the f(/3, K) 
by their average (1.6). The approximate evolution equation can be brought into the 
form: 

^f(l3, = (dV + m)m Ni) + VY, Cj( N i, V)f(/3, N ± + j). (2.1) 
In the combined limit (1.7) the coefficients Cj(Ni, V) were found to be 



(eqs(5.5-ll) of paper 1): 

c^ 2 (x) = - Wl , (2.2) 

c_i(x) = Sw! - w 2 - 2(2d - l)x, (2.3) 

co(x)= -3w 1 + w 2 + 4(2d - l)x - d + x, (2.4) 

Cl (x)= wi-2(2d-l)x + d-x. (2.5) 



Here wi(x) and w 2 (x) are the probabilities to find on the valence bond configurations 
K(Ni) dimer pairs with arbitrary and parallel orientation, respectively. The approxi- 
mate evolution equation (2.1) can be solved analytically in the combined limit (1.7) if 
we assume that the averages (1.6) have the form (1.8). Then eq.(2.1) yields a partial 
differential equation for <f>(/3,x): 

-AogR (2.6) 

L(x,R) (2.7) 
where 

L(x,R) = d + 3x + J2cj(x)R j (2.8) 

j 



d_ 

83 
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with the initial condition: 

(f)((3 = 0, x) = 0. (2.9) 

The solution proceeds as follows. We first introduce R instead of /3 — /3(x, R) as a 
new variable by means of the Legendre transform: 

$(x,i?)=#,x)-^,x). (2.10) 



Then we find for the partial derivatives: 

R — f: 

dx 



<9$ 

dx 



hgR- (3-^-L(x,R), (2.11) 



R 



<9$ 
OR 

If we eliminate (3: 



-pJL L (x,R). (2.12) 



3RaI = {^- logR )m (213) 

we get a partial differential equation for $(x, R) with initial condition: 

$(x,R = l) = 0. (2.14) 
A solution of eqs.(2.13), (2.14) is easily found: 

<S>(x,R)=x\ogR- J^-g(R',L(x,R)) (2.15) 
i 

provided that g(R',L(x,R)) solves the implicit equation: 

x = g{R' = R,L(x,R)). (2.16) 
Therefore we get the integrand g(R,L) in eq.(2.15) from 

L(g,R)=L. (2.17) 

3. Thermodynamical Properties in One Dimension 

The probability w±(x) to find dimer pairs on the valence bond configurations K(N 1 ) 
can be calculated analytically in the case d — 1, as is done in appendix A: 

w 1 (x) = (3.1) 

The implicit equation (2.17) for g(R,L) turns out to be quadratic and can be easily 
solved: 

g ± (R, L) — 1 — ^ (l - b(R) ± y/D(R, L)) (3.2) 

where 

a(R) =AR-11 + 5R~ 1 -R- 2 , (3.3) 
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b(R) = -4R+ U-8R- 1 + 2R- 2 , (3.4) 

D(R, L) — (L — b(R)) 2 - Aa{R) { -?-^-. (3.5) 

The solution has to satisfy eq.(2.16) which means in particular for: 

x = 0, L(x = 0, R) = R, g{R, L = R) = 0. (3.6) 

We have plotted in Fig.l both solutions given by eq.(3.2). The condition (3.6) is 
satisfied for: 

g + (R, L = R) = if 1< L < L x = 2.387425 (3.7) 

g-(R, L = R) = if L l <L<L 2 = 2.588229. (3.8) 

Both solutions meet each other at R — R^L), which defines the zero of D(R, L). This 
zero produces a branch point singularity of g(R, L) in the complex i?-plane, which has 
to be taken into account if we define the path of integration for the solution (2.15): 

<Z>(x,R) =x log R-J —g + (R',L) for L = L(x,R)<L l (3.9) 
i 



R2(L) R 

*(x,R)= xlogR- j —g + {R',L)~ j —g-{R',L) (3.10) 

1 ' ' R2{L) 

for Li < L(x, R) < L 2 

Finally, we need the relation (2.12) for going back to the original variable, the inverse 
temperature (5: 

/3{x,R) = J —— g+ {K,L) for L = L(x,R)<L 1 (3.11) 
i 



R2(L) R 

r HR' r> r HR' rt 

^• fl >= / 1FTl^ r '- l>+ I -WW-^ (312) 

1 R2{L) 

for Li < L(x, R) < L 2 . 

Since R and L vary only in open and finite intervals: 

1 < R < R 2 {L 2 ) = 3.772126, 1 < L < L 2 (3.13) 

the derivatives (2.6) and (2.7) are finite and the premise for eq.(1.9) is satisfied. There- 
fore the thermodynamical properties are obtained from the solution (f>(/3, x) at vanish- 
ing dimer density: 

x = 0, L = R= Jr0(/3, x = 0) = U(P). (3.14) 
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The inverse temperature /3(U) as function of the internal energy U can be taken from 
eqs. (3.11,12) for (x = 0, L = R = U). This function is shown in Fig. 2. It develops a 
singularity for 

R = U^L 2 , (3^oo. (3.15) 

which can be identified with the low temperature limit. Moreover we see from eq.(1.9) 
that Z/2 is just the ground state energy per site ho — L% — 2.588229. This value is 7% 
below the exact value h Q = 4 log 2 which gives us an idea on the accuracy of the lowest 
order in the dimer cluster expansion. In Fig. 3 we present the specific heat as function 
of the temperature^. Position and height of the maximum are: 

T = 0.982(0.962) C = 0.296(0.350) (3.16) 

For comparison we have listed in brackets the estimates of ref . [2] , obtained from small 
rings up to 16 sites. The low temperature behavior [3] of the specific heat, as it is 
predicted by conformal invariance and derived rigorously by J. Suzuki, Y. Akutsu and 
M. Wadati: 

C(T) = \t for T^0 (3.17) 

is represented in Fig. 3 by the dotted line. The high temperature expansion has been 
computed to order 21 [4] and is represented in Fig.3 by the dashed line. 



4. Thermodynamical Properties in Two and Three Dimen- 
sions 

The probabilities Wi(x),W2(x) to find on the valence bond configurations K(N\) dimer 
pairs with arbitrary and parallel orientation are determined in appendix A: 

Wl (x) = (2d-1) 2 -^- (4.1) 
a — x 

2 

W2 ( x ) = 2(d-l)j^-. (4.2) 

The implicit equation (2.17)-together with (2.8) and (2.2-5)- turns out to be of third 
order now: 

g 3 a(R, d) + g 2 (b(R, d) - 2-^) + g(c(R, d) + 3^) = L - dR (4.3) 

Ct (Jj 

where 

(2d-l) 2 1 1 
a(R, d) = 2 V R2(J2 J + —(-2Ad 2 + 18d - 4) + ^(24d 2 - lOd + 8) - 8R (4.4) 

{2d -I) 2 1 1 
b(R, d) = - D9 / + -j^(12d 2 - 2d - 1) - -(12d 2 + lOd + 5) + AR{d + 2) (4.5) 
R d dR d 

c(R,d) = 4(1 -2d) + 8d- R(4d + 2). (4.6) 
R 



J Our definition (1.1) and (1.3) of H and (3 differ from the usual one. Our (3 is related to the usual 
definition of the temperature via: (3 = (2dT)~ 1 . 
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The real solutions g±{R, L) are shown in Figs. 4 and 5 for d = 2 and d — 3, respectively. 
We observe that the condition (2.16) for x — 0, L(x = 0, i?) = <iR: 

9<7 (R,L = dR) = 0, a = +,- (4.7) 

is satisfied for the solution: 

a = + : if 1 < L < Li(d) (4.8) 

a = - : if Li(d) < L < L 2 (d) (4.9) 

where 

Li(d= 2) = 3.625833, L 2 (d = 2) = 4.154262 (4.10) 

Li(d= 3) = 4.862986, L 2 (d = 3) = 5.399962. (4.11) 

Both solutions meet each other at R — R2(L), where the derivative ^ divergies. At 
this point g(R, L) has a branch point singularity in the complex i?-plane which we pick 
up in the path of integration for the solution (2.15). The thermodynamical properties 
of the AFH model in d dimensions are obtained from: 

/" dR' d 

(3(x = 0, R) = J —— g+ (R>, L) for L = dR < L x (4.12) 
i 



R 2 (dR) 
- dR' d 



/a H a 
— —g + (#,L) (4.13) 

i 

+ J -^TqZ9-{R\L) for L x <dR<L 2 . 



dR' d 

R 2 (dR) 

This function relates the inverse temperature (3 with the internal energy per site: 



U([3) = ^0(/3, x = 0) = L = dR (4.14) 

which is shown in Fig. 6 for d — 2 and d = 3. A singularity appears for: 

dR->L 2 , [3(0,R)^oo (4.15) 

i.e. in the zero temperature limit. L 2 is the groundstate energy per site in the lowest 
order of the dimer cluster expansion (cf.(2.1)). The values listed in eqs. (4.10,11) have 
to be compared with the estimates of ref. [5,6,7] for the groundstate energy per site 
e(d): 

h (d = 2) = 2 - Ae(d = 2) = 4.678, h (d = 3) = 3 - Ae(d = 3) = 6.611 (4.16) 

The specific heats for d = 2 and d — 3 are shown in Figs. 7 and 8, respectively. 

For d = 2 the low temperature behavior is known from chiral perturbation theory [8]: 

= ^C(|) T 2 for T _^ / 417 x 
7 



where 

c=1.68 (4.18) 

is the spin wave velocity [5]. 

The high temperature expansion has been computed up to order 10 for d — 2 [9] 
and d = 3 [10]. The high and low temperature limits are represented in Fig. 7 by 
the dashed and dotted curves. We observe good agreement in the high temperature 
regime, whereas the low temperature behavior (4.17) cannot be reproduced by the 
lowest order of the dimer cluster expansion. 

5. Summary and Perspectives 

In this paper we have started a first attempt to evaluate the evolution equation (1.5) 
for the expectation values (1.4) of the Boltzmann factor between valence bond states. 
To lowest order in the dimer cluster expansion the thermodynamical limit of the evo- 
lution equation (2.1) leads to a partial differential equation (2.6-7) which can be solved 
analytically up to integrations. As a result, we obtain an explicit representation (3.11- 
12, 4.12-13) of the inverse temperature (3 as function of the internal energy per site 
U for the antiferromagnetic Heisenberg model in dimensions d = 1,2,3. The predic- 
tions of the lowest order in the dimer cluster expansion for the specific heat are in 
agreement with the high temperature expansion but they are in disagreement with 
the low temperature expansion (3.17) and (4.17). This means, that the formation of 
dimer clusters-which has not been taken into account in the lowest order approxima- 
tion (2.1)- plays an important role for the correct description of the low temperature 
behavior. The specific heat has a maximum at temperatures T(d) which decreases 
with the dimension d. The height of the maxima is almost independent of d, the 
width shrinks with d. Therefore, we do not see any indications for a phasetransition 
in our results for the specific heat in three dimensions, as it was suggested from the 
hightemperature expansion of ref.[10]. 
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Appendix 



A Dimer Pairs on Valence Bond Configurations 

In this appendix we will compute the probability: 

™([1,2][3,4]) = vimy 1 £ 5([1,2][3,4] G K(Nx)) (Al) 

if(JVi) 

to find a dimer pair at the neighbouring sites 2 and 3 on a valence bond configuration 
K(Ni) with Ni dimers. v(Ni) denotes the total number of these configurations. 
Let us keep the Ax dimers fixed in a 'dimer configuration' D(Nx) and ask how many 
possibilities we have to connect the remainig V — 2Nx sites on the lattice by valence 
bonds. This number is approximately given by (V — 2Ni — 1)!! i.e. by the total 
number of valence bond configurations on the remaining sites. In a strict sense we 
should count only those configurations on the remaining sites which do not contain 
any further dimers. Indeed these configurations dominate the thermodynamical limit, 
since configurations with dimers have lower entropy. This argument tells us that the 
dimer probability (A.l) on the valence bond configuration K(Ni) is just the same as 
on the dimer configuration D(Ni): 

w([l,2}[3,4}) = p(N 1 ,V)- 1 £ *([l,2],[3,4]e£>(JV-x)) (A2) 

D(Ni) 

where: 

p(N 1 ,V)= £ (A3) 

D(Ni) 

is the number of possibilities to distribute Nx dimers on a lattice with V sites. 
For d — 1 and open boundary conditions one finds: 

where n\ = V — 2Ni is the number of sites which are not occupied with dimers. 
For d — 1 and periodic boundary conditions one finds: 

p(N x , V) = Po (Nx, V) + p (A a - 1, V - 2) - poiNx + 1, V - 2). (A5) 

The second and third term take into account the different interpretation of the valence 
bond [1, V}. It is a dimer for periodic but not for open boundary conditions. 
In the case d — 1 the probability (A. 3) to find a dimer pair turns out to be: 

which leads in the combined limit (1.7) to eq.(3.1). 

In the case d = 2 and d = 3 the dimer pair probability (A. 2) can be determined in a 
simulation of the monomer dimer system. We have done that on the lattices |: 

2 The simulations on the 200 3 lattice has been performed on the connection Machine CM5. The 
results are based upon a test version of the software where the emphasis was on providing functionality. 
This software release has not had the benefit of optimization or performance. 
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• 100 x 100, 200 x 200 



• 50 x 50 x 50, 200 x 200 x 200 

The results for the dimer pair probabilities wi(x),l = 1,2 together with the fits 
(4.1-2) are shown in Fig. 9 for d = 2 and in Fig. 10 for d = 3. The behavior of 
wi(x) for a; = N\/V — > and x — > 1/2 and that of w 2 (x) for x — > can be easily 
checked by combinatorical considerations. Note that the fit (4.1) for the dimer pair 
probability Wi(x) is perfect for d = 2 and d = 3. We therefore believe, that (4.1) 
represents indeed the exact formula for this quantity. On the other hand, the fit (4.2) 
for the parallel dimer pair probability W2(x) is not perfect. We therefore expect that 
the exact formula for this quantity is more complicated than the fit (4.2). 
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Figure Captions 

1. The real solutions (3.2) of the implicit equations (2.17) for d — 1 . The solutions 
coincide on the solid curve R = RziL) where the square root in eq. (3.2) vanishes. 

2. The inverse temperature f3(U) versus the internal energy per site U, as it follows 
from (3.11,12) for x = 0, U = L = R, d = 1. 

3. The specific heat of the one-dimensional AFH- model in lowest order of the dimer 
cluster expansion. The dotted curve represents the low temperature behavior 
(3.17), the dashed curve the high temperature expansion of ref.[4]. 

4. The real solutions g±(R; L) of the implicit equation (4.3) for d = 2. The solutions 
coincide on the solid curve R = R2(L). 

5. Same as Fig. 4 for d — 3. 

6. The inverse temperature f3(U) versus the internal energy per site U as it follows 
from eqs. (4.12,13) for x — 0, U = L = dR. Solid curve: d = 2, dotted curve: 
d = 3. 

7. The specific heat of the d = 2 AFH model in lowest order of the dimer cluster 
expansion. The dotted curve represents the low temperature behavior (4.17); 
the dashed curve the high temperature expansion of ref. [9]. 

8. Same as Fig. 7 for d — 3. 

9. The dimer pair probabilities for d = 2: solid curve: Wi(x), dotted curve: w 2 (x) 
10. Same as Fig. 9 for d — 3. 
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